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Abstract 



pi i' Wc study a Hamiltonian quantum formalism of a spherically symmetric space-time which can 

D ' be identified with the interior of a Schwarzschild black hole. The phase space of this model is 

' spanned by two dynamical variables and their conjugate momenta. It is shown that the classical 

Lagrangian of the model gives rise the interior metric of a Schwarzschild black hole. We also show 
that the the mass of such a system is a Dirac observable and then by quantization of the model 
by Wheeler-DeWitt approach and constructing suitable wave packets we get the mass spectrum 
CJ ' of the black hole. 
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^ '. 1 Introduction 
m 

ff^ . Black hole physics has played a central role in conceptual discussion of general relativity in classical 

and quantum levels. For example regarding event horizons, space-time singularities and also studying 
QO ' the aspects of quantum field theory in curved space-time. In classical point of view, the horizon of a 

black hole which is a one way membrane, and the space-time singularities are some interesting features 
of the black hole solutions in general relativity [1]. In spirit of the Ehrenfest principle, any classical 
adiabatic invariant corresponds to a quantum entity with discrete spectrum, Bekenstein conjectured 
that the horizon area of a non extremal quantum black hole should have a discrete eigenvalue spectrum 
p]. Also, the black hole thermodynamics is based on applying quantum field theory to the curved 
\ space-time of a black hole [3]. According to this formalism, the Hawking radiation of a black hole 

is due to random processes in the quantum fields near the horizon. The mechanism of this thermal 
radiation can be explained in terms of pair creation in the gravitational potential well of the black 
hole d]. The conclusions of the above works are that the temperature of a black hole is proportional 
to the surface gravity and that the area of its event horizon plays the role of its entropy. In this 
scenario, the black hole is akin to a thermodynamical system obeying the usual thermodynamic laws, 
often called the laws of black hole mechanics, first formulated by Hawking [3]. In more recent times, 
this issue has been at the center of concerted efforts to describe and make clear various aspects of 
the problem that still remain unclear, for a review see [5]. With the birth of string theory [6], as 
a candidate for quantum gravity and loop quantum gravity [7], a new window was opened to the 
problem of black hole radiation. This was because the nature of black hole radiation is such that 
quantum gravity effects cannot be neglected [8]. According to all of the above remarkable works, it is 
believed that a black hole is a quantum mechanical system and thus like any other quantum system 
its physical states can be described by a wavefunction. Indeed, due to its fundamental conceptual 
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role in quantum general relativity, we may use it as a starting point for testing different constructions 
of quantum gravity [9l [10] . 

On the other hand, studying the interior of the Schwarzschild black hole is interesting for various 
reasons Many authors expressed the idea that the interior of a black hole can be considered as 
an anisotropic cosmological model [12]. In this direction, the author of [13] used it as a cosmological 
model to describe time dependent cosmological constant. One of the most important and interesting 
cosmological usage of the interior solution is considering black hole as a spawn of mother universe 
so that our universe was born from it as a daughter universe [H]. Many cosmological models in 
which our universe emerges from the interior of a black hole have been proposed in [15]. From a 
mathematical point of view, the space-time metric of the interior of a black hole can be constructed 
from its exterior metric. Indeed, one of the interesting features of general relativity is to generate 
cosmological solutions from known static solutions of the Einstein field equations. The methods from 
which one finds such solutions are investigated in |16j . Generally, these methods are based on a 
diffeomorphism between the known static solutions and the corresponding cosmological model. For 
example, the coordinate transformation t -H- r converts the Schwarzschild metric to a dynamical 
model which can be identified with its inside space-time. Indeed, at the horizon of a Schwarzschild 
black hole the light cone tips over and for r < Vg, d/dr become time-like while d/dt become space- 
like. This means that inside the black hole the metric has time-dependent coefficients. Thus, one can 
use this correspondence to make a quantum theory of black holes based on a quantized cosmological 
model. 

In this letter we deal with the Hamiltonian formalism of a time-dependent spherically symmetric 
space-time. We show that the classical solutions of such a dynamical system can be identified with 
the interior space-time of a Schwarzschild black hole. In this model the metric functions play the role 
of independent dynamical variables which with their conjugate momenta construct the corresponding 
phase space. We then consider a Hamiltonian quantum theory by replacing the classical phase 
space variables by their Hermitian operators. We study the various aspects of the resulting quantum 
model and corresponding Wheeler-DeWitt (WD) equation and present closed form expressions for the 
wavefunction of the black hole. We shall see that the quantum solutions also represent a quantization 
rule for the mass of the black hole. In what follows, we work in units where c = h = 1. 



2 The model 



We start with the Einstein-Hilbert action 
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(fxy/^TZ + Sygh-, 



(1) 



where TZ is the Ricci scalar, g is the determinant of the metric tensor and Sygh is the York-Gibbons- 
Hawking boundary term . We assume that the geometry of space-time is described by a spherically 
symmetric metric with time-dependent line element as 



N\t) 



dt"^ + u{t)dr^ + h^{t){d'd^ + sin^ M^"^ 



(2) 



where N{t) is the lapse function, while z^(t) and h{t) are functions of t only, which play the role of our 
dynamical variables to construct the phase space. It is clear that the above metric with t = const, and 
r = const, describes the surface of a 2-sphere with radius h{t) and area A = 47r/i^(t). By substituting 
([2|) into ([T|) and integrating over spatial dimensions, we are led to an effective Lagrangian in the 
minisuperspace (z/, h) as 



Vo 
'SttG 



(3) 
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in which Vq is the volume of the space part of the action where is treated to be a finite constant. A 
point {u, h) in this minisuperspace represents a 4-geometry. Although, in this step we can vary the 
above Lagrangian to get the equation of motion for h and but the Hamiltonian constraint resulting 
from this Lagrangian does not have the desired form for construction of WD equation describing the 
relevant quantum model. Thus, to transform Lagrangian ([3]) to a more manageable form, consider 
the following change of variables 

X — y = h, X -\- y = hv. (4) 
In terms of these new variables, Lagrangian ([3]) takes the form 



(5) 



SttG 

Now, if we use the following coordinate transformation which is introduced in 

(x - y)i = i(u - v), 

(6) 

{x + y)2 = l(u + v), 

and the lapse rescaling 

N = {x^ -y'^)^N, (7) 

then the corresponding Lagrangian becomes 

C = -MliV^N-^ii? - v^) + M^iVoN{u^ - v^), (8) 

which describes an isotropic oscillator-ghost-oscillator system. 

To construct the Hamiltonian of the model, note that the momenta conjugate to u and v are 

_ 2yoM| 



(9) 

Also, the primary constraint is given by 



% = ^ = 0. (10) 
dN 

In terms of the conjugate momenta the Hamiltonian is given by 

H = iiUu + vlly - C, (11) 

leading to 



H = -7T7JTF(n' - n^) - VoM],^Niu' - v^). (12) 



Because of the existence of constraint (jlOp . the Lagrangian of the system is singular and the total 
Hamiltonian can be constructed by adding to H the primary constraints multiplied by arbitrary 
functions of time X{t) 

Ht = -^]^^(n« - n?,) - VoMl^iNiu^ - v^) + XU^. (13) 
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The requirement that the primary constraint should hold during the evolution of the system means 
that 



UN = {nN,HT}-0, (14) 
which leads to the secondary (Hamiltonian) constraint 

^ = -J^T^jri^l - n^) - ^oM2,(n2 - v^) « 0. (15) 

Before going any further, some remarks are in order about the mass of the system. In the usual 
classical theory, the unique solution to the vacuum Einstein equation for the spherically symmetric 
space-time is the Schwarzschild solution in which there exists an integration constant representing 
the mass of the black hole. In the canonical formalism for the spherically symmetric space-time, by 
following Fischler-Morgan-Polchinski |17j and Kuchar [18] calculations the mass of a black hole can 
be regarded as a dynamical variable, which can be expressed as a function of canonical data. The 
spherically symmetric hypersurface on which the canonical data is given is supposed to be embedded 
in a Schwarzschild black hole space-time whose metric is given by 

2MG , , , 

N = \, ^ = — 1. h = t. (16) 

This identification of the space-time with the canonical data enables us to connect the Schwarzschild 
mass M with the canonical data on any small piece of a space- like hypersurface. As a result, from 
combination of the relations dH), dS]), (l9|) and (I16p the mass function can be expressed by the canonical 
variables as 

= {j^, (n. + nj' + v„Mh(u - vf) . (17) 

It is easy to show that the Poisson bracket of mass function with the Hamiltonian vanishes strongly 

{M,n} = 0, (18) 
which shows that the mass is a constant of motion. 



3 Classical solutions 

The setup for constructing the phase space and writing the Lagrangian and Hamiltonian of the model 
is now complete. The classical dynamics is governed by the variation of Lagrangian ([8]) with respect 
to u and v, that is 

u + N'^u = 0, v + N'^v = 0. (19) 

Up to this point the model, in view of the concerning issue of time, has been of course under- 
determined. Before trying to solve these equations we must decide on a choice of time in the theory. 
The under-determinacy problem at the classical level may be removed by using the gauge freedom 
via fixing the gauge. A glance at the above equations shows that a suitable gauge is = Mpi which 
results the following equations 

il + Mliu = 0, v + Mliv = 0. (20) 
Choosing some parameters 9i and 62 as the integration constants, the solutions are obtained as 

u = Acos{Mpit + ei), V = Bcos{Mpit + e2), (21) 
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where A and B are constants to be determined later. Now, the above solutions should satisfy the 
constraint of a vanishing Hamiltonian. Thus, substitution (j2ip into (jl5p gives a relation between the 
constants A and B 

A = ±B, (22) 



implying that we can rewrite the solutions (I2ip as 

u = Acos{Mpit + 0i), V = Ar]cos{Mpit + e2), (23) 

where rj takes the values ±1 according to the choices in (I22p . Since in the quantum version of the 
model we are interested in constructing wave packets from the WD equation, we would like to obtain 
a classical trajectory in configuration space {u,v), where the classical time t is eliminated. This is 
because no such parameter exists in the WD equation. It is easy to see that the classical solutions 
(j23p may be displayed as the following trajectories 

u'^ + v^ - 2ri cos{M)uv - A^ sin^{M) = 0, (24) 

where = 9i — 62- Equation (I24|) describes ellipses which their major axes make angle 7r/4 with 
the positive/negative u axis according to the choices ±1 for rj. Also, the eccentricity and the size of 
each trajectory are determined by and A respectively. 

On the other hand, inserting solution (j2ip into the mass function, we obtain 

M = 2'kMIiA^ sin(— ). (25) 

To obtain the usual Schwarzschild solution, consider time rescaling Ndt = dr, where from d?]) we 
have A'" = Mpi{u^ — u^)/4. Consequently, solutions (pTj) give 

cos(2Mp;f + ^1 + 62) = ~ .^"^ , (26) 

A'^ sm l\U 

where tq is a constant of integration. From equations and ([6]) one then has 



-2coth(4^)r, ?? = 1, A6'<0, 
2tanh(4^)r, r? = -1, A6I > 0. 



h={.:Tie\"' " : :::: (27) 



and 



'cot2(M)(2CiM,^t(M)_i), ^ = 1, A^<0, 

tan2(M) ( 2^ tan(^) - l) , r? = -1, A^ > 0, ^ ' 



in which we put 



To = -lA^s\nM. (29) 

Thus, in terms of the metric functions (/i, u) the classical trajectories ()24p take the form 

_l / A^siii^ M \l n + rjcosAeX 
^ ~ 2 Vl -r?cosA6' ) h~ Vl -??cos A^ ) ' 

Therefore, the above Lagrangian formalism leads us to the following interior Schwarzschild black 
hole metric if we assume A^ = 2n7r — 7r/2 when ij = 1 and A6 = 2n7r + 7r/2 when r] = —1, respectively 

ds^ = _ _ 1^ dr' + _ 1^ dr^ + t2 {d^^ + sin2 Mip') . (31) 
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It is clear from the condition i/(r) > that this solution is only valid for r < 2MG. The above 
metric has an apparent singularity at r = 2MG. This singularity is like the coordinate singularity 
associated with horizon in the Schwarzschild or de Sitter space-time, and as is well known, there 
are other coordinate system for which this type of singularity is removed [1]. Another singularity 
associated with the metric (13ip is its essential singularity at r = 0. As we know, in general relativity, 
to investigate the types of singularities one has to study the invariants characteristics of space time 
and to find where these invariants become infinite so that the classical description of space-time 
breaks down. In a 4- dimensional Riemannian space-time there are 14 independent invariants, but to 
detecting the singularities it is sufficient to study only three of them, the Ricci scalar TZ, R^yR^'^ and 
the so-called Kretschmann scalar R^^^jsR^'^'^^ . For the metric ()3ip the Kretschmann scalar reads 

K = R,,,sR^-'^' = ^^^. (32) 

Now, it is clear that the space-time describing by the metric (j3ip has an essential singularity at r = 0, 
where can not be removed by a coordinate transformation. Note that the Schwarzschild manifold 
contains an anisotropic expanding universe, the white hole portion of the extended geometry, and also 
an anisotropic collapsing universe, the black hole interior as well. In this paper we focus attention on 
the black hole interior portion of the geometry, but all conclusions may be restated in terms of the 
expanding white hole geometry due to the time reversal symmetry of both Schwarzschild geometries. 



4 Dirac Observables 

As is well known general relativity is invariant under the group of diffeomorphisms of the space-time 
manifold M. The main consequences of such a diffeomorphism invariance are that the Hamiltonian 
can be expressed as a sum of constraints and any observable must commute with these constraints. 
An observable is a function on the constraint surface such that is invariant under the gauge trans- 
formations generated by all of the first class constraints. By a first class constraint we mean a phase 
space function with the property that it has weakly vanishing Poisson bracket with all constraints. 
As an example the momentum and Hamiltonian constraints are always first class, see (1101) and (jlSp . 
The Hamiltonian and momentum constraints in general relativity are generators of the corresponding 
gauge transformations, and so a function on the phase space is an observable if has weakly vanishing 
Poisson brackets with the first class constraints. To find gauge invariant observables, we can proceed 
as follows. In Lagrangian dH]), the unconstrained phase space P is with global canonical coordinates 
(xi,nj), i = 1,2, with Poisson brackets {xi^Hj} = 5ij. Now, let us define on P the complex-valued 
functions 

(33) 

c* = MpiVVox. - i^j^n,, 

where Xj = {u,v} and Hj = {n„,n„}. The set S = {Cj,C*,l} on P is closed under the Poisson 
bracket, {Cj,C|} = —iSij and every sufficiently regular function on P can be expressed in terms of 
the sums and products of the elements of 5. Hence, the Hamiltonian and mass can be viewed as 

7i = C*C^ — C*Cu, (34) 

and 

M = ^^{Cu-C:){C:-C,). (35) 

Therefore, the classical Poisson algebra generated by the elements of 5 is sufficiently large for describ- 
ing the classical dynamics of the system. In the next section we will use this algebra as the starting 
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point for quantization of the model. The classical dynamics of these variables is 

(36) 

To find a set of constraints of motion, consider on T the fmictions 

Jq = —^{C*Cv + C*Cu), 

J_i_ = CyCu, (37) 

■J- — '^V^UI 

whose Poisson brackets form a closed algebra by 

{J+,J_} = 2iJo, 

(38) 

{Jo, J±} = TiJ±- 

The elements of this algebra have strongly vanishing Poison brackets with the Hamiltonian. Since 
the physical space is two dimensional, there will be at most two independent constraints. On the 
constraint surface Ti = 0, the functions J's are not algebraically independent but satisfy the identity 

4 - J+J- = 0. (39) 
Also, we can write the mass function as a combination of J's as 

M = -^(2Jo + J+ + J_). (40) 
Thus, the quantity M is a gauge invariant observable. 

5 Quantization of the model 

We now focus attention on the study of the quantization of the model described above. The quan- 
tum version of the model described by relations {xi,Ilj} = 5ij can be achieved via the canonical 
quantization procedure which leads the following commutation relations 

[u,Uu]=i, [v,U^] = i. (41) 

Then, the set of hermitian quantum operators S = {Ci,Cj, 1} will have the following commutator 
algebra 

[Ci,C]]=6i,l, [C,,1] = [C/,1]=0. (42) 

The set S and its commutator algebra are the quantum counterpart of the set S and its Poisson 
bracket algebra. The operator version of the classical Hamiltonian constraint takes the form 

H = ClC - ClCu. (43) 

Let us define vacuum state according to 

C„|0,0>=0, a|0,0>=0, (44) 
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so that Cj and Cj are annihilation and creation operators respectively. Note that by the Hamiltonian 
operator (03]), zero point energy is canceled out. Consequently, the WD equation can be written as 



(45) 



where oj = Vq/8ttG. This equation is a quantum isotropic oscillator-ghost-oscillator system with 
zero energy. Therefore, its solutions belong to a subspace of the Hilbert space spanned by separable 
eigenfunctions of a two-dimensional isotropic simple harmonic oscillator Hamiltonian. Separating the 
eigenfunctions of (j45p in the form 



(46) 



yields 



^1/4 
,V4 



V2"ini! 



-uju'^/2 



(47) 



subject to the restriction ni = 1x2 = n. In (|47p . Hn[x) are Hermite polynomials and the eigenfunctions 
are normalized according to 



/+00 
e-""" Hn{x)Hm{x)dx = 2"7ri/2n!J^,. 
-00 



(48) 



The set {$ 

ni,n2(^'^)} forms a closed span of the zero sector subspace of the Hilbert space L2 of 
measurable square-integrable functions on with the usual inner product defined as 



/ 



^ni,n2 {u, 'V)$mi,m.2 {u, v)dudv 



n2,m2 ) 



(49) 



that is, the orthonormality and completeness of the basis functions follow from those of the Hermite 
polynomials. Therefore, in the position representation, we may write the general solution of the WD 
equation as a superposition of the above eigenfunctions 



n=0 



(50) 



where Cn are a set of complex constants. In general, one of the most important features in quantum 
WD approach is the recovery of classical solutions from the corresponding quantum model or, in 
other words, how can the WD wavefunctions predict a classical model. In this approach, one usually 
constructs a coherent wave packet with good asymptotic behavior in the minisuperspace, peaking in 
the vicinity of the classical trajectory. Therefore, for our subsequent analysis, by using the equality 



V —Hn{x)Hn{y) = exp 

n=0 



2xys — s (X + y 
2(1 -s2) 



(51) 



we can evaluate the sum over n in (j50p and simple analytical expression for this sum is found if we 
choose the coefficients Cn to be c„, = i?2" tanh^, where B and ^ are arbitrary complex constants, 
which results in 



*(x, y) = 7Vexp[-| cos(2/3) cosh(2a)(a;2 + y^ - 277tanh(2Q)xy)] 
X exp[- J sinh(2a) sin(2/3)(x^ + 2/^ - 2r? coth(2a)xy)]. 



(52) 
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Figure 1: The figure shows the square of the wavefunction. This figure is plotted for 

numerical values a = 1, 6a = 0.3, /3 = 7r/4, Sf] = 37r/20, r] = 1 and u = 1. 

In this expression a and f3 are the real and imaginary parts of C respectively, that is, C = a + 
X = ^/uiu, y = ^/ujv and M is a numerical factor. In figure [T] we have plotted the square of the 
wavefunction for typical values of the parameters in which we have taken the following combination 
of the solutions 

<^{x, y) = ^r,,aA^^ y) - ^■/y,a+<5a,/3+<5/3, (53) 

for some 5a and 5/3 in the vicinity of a and (3. It is seen that a good correlation exists between the 
quantum pattern shown in this figure and the classical trajectories (j24p in configuration space (u, v). 

In continuation of this section we focus attention on the Dirac observables of the black hole and 
try to build an algebra of physical operators, that is, a subalgebra of A that would leave Vhh,o 
invariant. To do this, notice that the Poisson bracket algebra of the classical J's can be promoted 
into a commutator algebra version by setting 

J+ = C,Cu, J- = Clcl J^ = -]^(clC, + ClCu + l) , (54) 
so that the corresponding commutators are 

[J+,J_] = 2Jo, [Jo,J±]=tJ±. (55) 

Equations (|55p are recognized as the commutators of the Lie algebra of 50(2,1). Since all these 
operators commute with Hamiltonian H, we choose our physical operator algebra to be the algebra 
generated by the set {J±,Jo,l}. Note that in analogy to the classical algebraic identity ([39|) the 
quantum J's are not independent and satisfy the following algebraic relation 

Jl-\{J^J^ + J-J+) + \ = ^. (56) 

Now, the action of the operators {Cj, C|} on the states of the physical Hilbert space is 

Cn|ni,n2 >= ^/nl\nl - l,n2 >, Ci\ni,n2 >= Vf^i + M^i + l,?i2 >, 

(57) 

Cv\ni,n2 >= ^/n^\nl,n2 - 1 >, Cj|ni,n2 >= ^/^^2 + l\ni,n2 + 1 >, 
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so that with the help of (|54|) . one can put the action of J's on the physical states as 



Jo|ni, n2 >= -2(^1 + "-2 + l)!"-!, "-2 >, 

J+\ni,n2 >= ^/nrn2\nl - l,n2 - 1 >, (58) 



J-\ni,n2 >= (ni + l)(n2 + l)|ni + l,n2 + 1 >, 

with constraint ?ii = 112- Consequently the expectation value of mass operator for the interior solution 
win be0 



which means that black hole is quantized in discrete states with energies E^, that is (in ordinary 
units) 

It is well known that the parameter M in the outside of a black hole is measured at spatial infinity 
where the Newtonian approximation is valid and so it is a measure of mass of the black hole. In this 
sense, we know that the energy spectrum of a quantized black hole from an outside observer point of 
view is (see [2] and pO] ) 

En = Mpic^V^, (61) 
which gives the following relation for the difference between two nearby states 

AE= JfjL'^ (62) 

\/n + 1 + vn 

Now, it is obvious that AE tends to zero as n — t- 00 which is in agreement with the correspondence 
principle and the black hole can be described classically in this regime. On the other hand, the physical 
meaning of black hole mass parameter may be different inside the event horizon. As equation (j60p 
shows, the energy eigenvalues consist of equally spaced spectrum which resemble the energy spectrum 
of the harmonic oscillator. This is not surprising, since a decomposition of the quantum fields inside 
the black hole into normal modes is essentially a decomposition into harmonic oscillators that are 
decoupled. Also, from (j60p we see even the lowest energy level, i.e. the level n = 0, has some nonzero 
energy, the so called ground state energy. The existence of such a energy level is a purely quantum 
mechanical effect and may be interpreted in terms of the uncertainty principle. Indeed, it is this zero- 
point energy that is responsible for the fact that the system under consideration (interior of the black 
hole in our case) does not "freeze" at extremely low energies. Another feature of the result (j60p is its 
agreement with the correspondence principle in the sense that — g-*-^ — >• 1 as n — )• 00. As the quantum 
number n increases the system tends to its classical regime which in this case is a superposition of 
the eigenstates of the type (|50p . Such a so-called coherent state consists of unblurred wave packets 
which minimize the uncertainty relation. 

6 Conclusions 

It is interesting to investigate how the interior of a black holes would be quantized. Discrete spectra 
arise in quantum mechanics in the presence of a periodicity in the classical system, which in turn leads 



^Note that according to the definition of volume Vo = 4% J dr, its dimension is length. 
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to the existence of an adiabatic invariant or action variable. Boher-Somerfeld quantization implies that 
this adiabatic invariant has an equally spaced spectrum in the semi-classical limit. Using this approach 
one can determine the mass and the area spectrum of the black holes in view of an outside observer [19] . 
In this paper we have evaluated analytically the mass spectrum of the interior Schwarzschild black 
hole by implementing WD equation. The corresponding phase space of the interior of a Schwarzschild 
black hole is spanned by two dynamical variables and their conjugate momenta. We have shown that 
the classical Lagrangian of the model gives rise the interior of Schwarzschild solution. Then, we 
studied quantization of the model through WD equation and by imposing suitable conditions on the 
solutions of WD equation we have obtained the mass spectrum of the interior of the black hole. Our 
result shows that the mass spectrum, not only is discrete, but also is equally spaced. 
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